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The  pulsed-field-ionization  zero-kinetic-energy  photoelectron  spectrum  of  Xe2  has  been  measured 
between  97  350  and  108  200  cm”1,  following  resonant  two-photon  excitation  via  selected 
vibrational  levels  of  the  C  0*  Rydberg  state  of  Xe2.  Transitions  to  three  of  the  six  low-lying 
electronic  states  of  Xe2  could  be  observed.  Whereas  extensive  vibrational  progressions  were 
observed  for  the  transitions  to  the  1(3/2 g)  and  1(3/2 u)  states,  only  the  lowest  vibrational  levels  of 
the  11(1  Hu)  state  could  be  detected.  Assignments  of  the  vibrational  quantum  numbers  were  derived 
from  the  analysis  of  the  isotopic  shifts  and  from  the  modeling  of  the  potential  energy  curves. 
Adiabatic  ionization  energies,  dissociation  energies,  and  vibrational  constants  are  reported  for  the 
1(3  Hg)  and  the  1(3  Hu)  states.  Multireference  configurational  interaction  and  complete  active  space 
self-consistent  field  calculations  have  been  performed  to  investigate  the  dependence  of  the  spin-orbit 
coupling  constant  on  the  intemuclear  distance.  The  energies  of  vibrational  levels,  measured 
presently  and  in  a  previous  investigation  (Rupper  et  al,  J.  Chem.  Phys.  121,  8279  (2004)),  were 
used  to  determine  the  potential  energy  functions  of  the  six  low-lying  electronic  states  of  Xe2  using 
a  global  model  that  includes  the  long-range  interaction  and  treats,  for  the  first  time,  the  spin-orbit 
interaction  as  dependent  on  the  intemuclear  separation.  ©  2008  American  Institute  of  Physics. 
[DOI:  10.1063/1.2937133] 


I.  INTRODUCTION 

The  rare  gas  dimer  ions  represent  model  systems  to 
study  ion-atom  collisions.  They  are  a  source  of  absorption 
losses  in  excimer  and  ion  lasers1-3  and  are  encountered  in 
high-pressure  lamps  and  plasmas.4  The  potential  energy 
curves  of  rare  gas  dimer  ions  are  the  key  to  understanding 
their  decay  and  fragmentation  dynamics,5-11  and  are  required 
to  model  highly  electronically  excited  states  of  the  neutral 
dimers  by  multichannel  quantum  defect  theory.  Detailed 
knowledge  of  the  potential  energy  functions  of  these  ions  is 
also  important  for  a  fundamental  understanding  of  the 
chemical  binding  and  the  fragmentation  dynamics  of  larger 
rare  gas  cluster  ions,12-22  and  is  a  prerequisite  for  the  calcu¬ 
lation  of  charge-exchange23  and  differential  scattering  cross 
sections24  that  are  necessary  to  model  the  effects  of  rare  gas 
plasmas  associated  with  electric  space  propulsion  thrusters. 

The  homonuclear  rare  gas  dimer  ions  RgJ  (Rg=Ne,  Ar, 
Kr,  Xe,  and  Rn)  possess  six  low-lying,  closely  spaced  elec¬ 
tronic  states.  At  short  intemuclear  distances,  these  states  can 
be  labeled  using  Mulliken’s  notation25  as  A  22{/2m,  B  2nftg, 
C  2nllw  and  D  2S{/2^  (D=  1 /2,3/2).  At  large  intemuclear 
distances,  the  spin-orbit  interaction  mixes  2  and  II  electronic 
states  of  the  same  value  of  D  and  same  g/u  symmetry.  The 
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appropriate  labeling  then  is  1(1 /2m),  1(3 /2g),  1(3 /2m), 
1(1  I2g),  11(1 /2m),  and  H(l/2g),  where  the  value  of  11  and 
the  g/u  symmetry  are  given  in  parentheses  and  the  labels  1 
and  II  designate  states  correlating  with  the  first  Rg(]50) 
+  Rg+(2P3/2)  and  second  Rg(150)-»-Rg'b(2/>1/2)  dissociation 
limits,  respectively. 

Whereas  early  investigations  were  primarily  concerned 
with  the  qualitative  features  of  the  electronic  structure  of  the 
rare  gas  dimer  ions,25  recent  experimental  and  theoretical 
studies  have  provided  quantitative  information  on  the  poten¬ 
tial  energy  curves  of  several  electronic  states  of  Net,26,27 
A^,27-29  Kr2, 27,30  and  XeJ.31,32 

Most  experimental  data  on  the  low-lying  electronic 
states  of  XeJ  have  been  obtained  by  He  I,33,34  threshold35-38 
and  pulsed-field-ionization  zero-kinetic-energy.32,39  (PFI- 
ZEKE)  photoelectron  spectroscopy,  and  electron  impact 
ionization.11  Previous  experimental  studies  of  the  photoelec¬ 
tron  spectrum  of  Xe2  have  provided  data  on  the  1(1 /2m) 
<-X0+  I(3/2g)«— X 0^  and  11(1/2 m)<-X0+  photoionizing 
transitions.  The  spectroscopic  data  available  on  the 

1(3 /2m) <— X 0*,  I(l/2g)4-X0;  and  U(l/2g)4-X0;  transi¬ 
tions  are  much  less  extensive.  The  vibrational  assignment  of 
the  1(1 /2m),  I(3/2g),  and  I(3/2m)  states  are  unambiguous. 
Lu  et  a/.38  suggested  the  origin  of  the  11(1 /2m)  state  to  be 
located  at  107  109  cm”1  above  the  neutral  ground  state, 
while  Rupper  et  al.  found  its  origin  to  be  located  at 
107  157.7  cm”1.  Tonkyn  and  White39  and  Rupper  et  al ,32 
have  determined  three  vibrational  levels  (u+= 0,1,2)  of  the 
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1(3/2 u)  state  by  PFI-ZEKE  photoelectron  spectroscopy  us¬ 
ing  single-photon  excitation  from  the  X  0*  ground  neutral 
state  of  Xe2.  The  1(1  /2g)  state  has  not  been  observed  experi¬ 
mentally  so  far  and  is  reported  to  be  either  repulsive4041  or 

H  31  4?  tq 

weakly  bound.  ’  ’  ’  Lu  et  al.  attributed  a  line  at 
—  1 07  715  cm-1  in  their  threshold  photoelectron  spectrum  to 
the  origin  of  the  11(1 12g)  state.  Rupper  et  al ?2  also  observed 
an  isolated  line  at  108  132.3  cm-1  in  their  spectrum,  which 
they  tentatively  attributed  to  the  II(l/2g)  (i/=0)«—  X0* 
transition.  When  combined  with  the  results  of 
semiempirical25  and  ab  initio  quantum  chemical31’40-43  cal¬ 
culations,  these  experimental  results  provide  a  qualitative  un¬ 
derstanding  of  the  relative  role  of  short-range,  long-range, 
and  spin-orbit  interactions  and  a  qualitative  description  of 
the  strengths  and  lengths  of  the  bonds  in  all  six  low-lying 
electronic  states  of  XeJ.  However,  despite  the  considerable 
experimental  and  theoretical  efforts  that  have  been  invested 
to  derive  accurate  potential  energy  functions  for  these  states, 
several  assignments  remain  controversial.  Moreover,  all  the¬ 
oretical  and  experimental  determinations  of  the  potential  en¬ 
ergy  functions  of  the  six  lowest  electronic  states  of  Xe£  have 
been  based  so  far  on  the  assumption  that  the  spin-orbit  cou¬ 
pling  constant  is  independent  of  the  intemuclear  separation 
R.  This  assumption  was  first  made  and  justified  by  Cohen 
and  Schneider44  for  Ne^,  but  has  never  been  rigorously  tested 
for  the  heavier  rare  gas  dimers. 

The  main  objectives  of  the  present  work  were  (1)  to 
obtain  a  complete  set  of  experimental  data  on  the  vibronic 
energy  level  structure  of  Xe£,  (2)  to  remove  all  persisting 
ambiguities  in  the  spectral  assignments,  and  (3)  to  identify, 
and  if  possible  to  quantify,  the  effects  of  a  potential  depen¬ 
dence  of  the  spin-orbit  coupling  constant  on  the  intemuclear 
distance  by  combining  the  results  of  high-resolution  photo¬ 
electron  spectroscopy  with  those  of  ab  initio  quantum 
chemistry. 

To  this  end,  a  new  investigation  of  the  1(3 /2g),  1(3/2 m), 
and  11(1  Hu)  states  of  Xe£  by  resonance-enhanced  two- 
photon  ionization  via  selected  vibrational  levels  of  the  C  0* 
intermediate  state  of  Xe2  was  carried  out.  Xe£  was  selected 
because  the  very  large  spin-orbit  interaction  seemed  ideally 
suited  for  the  detection  of  a  potential  /^-dependence  of  the 
spin-orbit  coupling  constant.  This  experimental  study  was 
accompanied  by  an  ab  initio  quantum  chemical  study  of  the 
/^-dependence  of  the  spin-orbit  coupling  constant. 

A  (1  + 1')  two-photon  excitation  scheme  via  selected  vi¬ 
brational  levels  of  the  C  0*  Rydberg  state  of  Xe2  located  in 
the  vicinity  of  the  Xe*([5p]56s'[l /2]1)+Xe(150)  dissocia¬ 
tion  limit  has  been  used  to  access  the  ionic  levels.  This 
species-selective  excitation  scheme  enables  one  to  avoid  un¬ 
desirable  contributions  to  the  spectra  originating  from  the 
ionization  of  free  xenon  atoms  or  of  large  clusters  which 
completely  obscure  several  regions  of  the  single-photon  PFI- 
ZEKE  photoelectron  spectrum.32  Moreover,  selecting  differ¬ 
ent  vibrational  levels  of  the  CO*  intermediate  state  enables 
one  to  access  different  vibrational  levels  of  the  cation.  Tran¬ 
sitions  to  the  vibrational  levels  up  to  u+=52  of  the  1(3/2 g) 
state  and  up  to  u+  =  22  of  the  1(3/2 u)  state  have  been  ob¬ 
served  and  unambiguously  assigned.  The  assignment  of  the 
11(1 /2m)  state  proposed  by  Rupper  et  al.32  was  confirmed  in 


an  analysis  of  the  isotopic  shifts.  In  the  range  between 
97  800  and  97  900  cm-1,  several  features  could  be  observed 
in  the  spectra  which  may  be  members  of  the  1(1  /2g) 
<-X0j  progression.  However,  an  unambiguous  assignment 
was  not  possible.  No  transitions  to  the  11(1 12g)  state  were 
observed  following  resonant  two-photon  excitation  via  the 
C  0*  Rydberg  state. 

II.  COMPUTATIONAL  DETAILS 

The  nuclear  charge  of  xenon  (Z=54)  requires  scalar  and 
spin-orbit  relativistic  effects  to  be  taken  into  consideration  in 
all  quantum  chemical  calculations  on  XeJ.  Therefore,  the 
scalar-relativistic  second-order45  and  eighth-order  (optimum 
parametrization)46  Douglas-Kroll-Hess  one-electron 
Hamiltonians47-49  were  employed  in  all-electron  calcula¬ 
tions.  Spin-orbit  coupling  was  considered  as  a  perturbation. 
For  the  approximation  of  the  total  electronic  wave  function, 
several  configuration-interaction-type  wave  functions  were 
investigated. 

The  first  model  chosen  was  a  minimal  multiconfigura¬ 
tion  self-consistent  field  wave  function  as  provided  by  the 
complete  active  space  self-consistent  field  (CASSCF) 
method.50  An  active  space  of  6  orbitals  (corresponding  to  the 
5/7  shell  of  Xe)  with  1 1  electrons  was  chosen  to  represent  the 
lowest  gerade  and  ungerade  electronic  states  with  either 
£  or  n  symmetry.  The  dynamic  electron  correlation 
was  estimated  in  a  multiconfiguration  second-order 
perturbation-theory51  (CASPT2)  approach  involving  the  4 d, 
5s,  and  5 p  shell  of  Xe.  State-average  CASSCF,  state-average 
CASPT2  and  state-interaction  spin-orbit  calculations  were 
carried  out  with  the  MOLCAS  electronic  structure  package53 
using  the  arbitrary -order  DKH  module.54’55  In  all  CASSCF/ 
CASPT2  and  spin-orbit  state-interaction  calculations  the 
atomic  natural  orbital  relativistic  core  correlated  (ANO- 
RCC)  basis  set  of  Roos  et  al.  of  the  following  size 
(22s/ 19/7/ \3d/5f/3g)/[\0s/9p/%d/5fl3g]  was  employed, 
where  the  figures  in  parentheses  refer  to  the  exponents  of  the 
basis  set  and  those  in  brackets  represent  the  contracted  basis 
functions.56  Atomic  mean-field  integrals57,58  (AMFIs)  and 
the  electronic  energy  obtained  from  the  CASPT2  calcula¬ 
tions  were  used  to  set  up  the  spin-orbit  matrix  for  the  state- 
interaction  procedure  in  order  to  obtain  the  spin-orbit 
coupled  states. 

To  assess  the  accuracy  of  the  minimal  CASSCF  model, 
we  performed  second-order  Douglas-Kroll-Hess  multirefer¬ 
ence  configuration  interaction  calculations  with  single  and 
double  excitations59  (MRCISD)  with  the  MOLPRO  program 
suite60  and  the  arbitrary-order  DKH  module,54  explicitly  cor¬ 
relating  35  electrons  (i.e.,  the  4 d,  5s,  and  5/7  shell  of  Xe).  An 
estimate  of  the  contribution  of  higher-order  excitations  to  the 
electronic  energy  was  made  by  taking  the  Davidson  correc¬ 
tion  into  account.61  In  the  MRCI  calculations  we  also  em¬ 
ployed  the  ANO-RCC  basis  set  of  Roos  et  al.,56  but  in  a  fully 
decontracted  manner.  All  calculations  were  carried  out  for 
intemuclear  distances  ranging  from  4.0  to  14.0  a.u.  In  order 
to  estimate  the  effect  of  core  electron  correlation  on  spectro¬ 
scopic  parameters  we  additionally  performed  Davidson- 
corrected  MRCISD  calculations  correlating  solely  the  va- 
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TABLE  I.  Computational  details,  equilibrium  bond  distances  Re,  and  harmonic  vibrational  constants  u>t  of  the  22*  and  1(1  Hu)  states.  The  all-electron  MRCI 
and  CASPT2  data  were  obtained  with  a  second-order  scalar-relativistic  Douglas-Kroll-Hess  (DKH2)  Hamiltonian.  The  symbol  Q  denotes  Davidson-corrected 
MRCI  results.  The  abbreviation  SOC  indicates  the  type  of  integrals  used  to  construct  the  spin-orbit  matrix  for  the  perturbative  state-interaction  approach 
(indicated  as  “+SO”  in  the  denomination  of  the  method).  For  further  information  on  the  methods  see  text.  The  coupled-cluster  CCSD-T  data  taken  from  Ref. 
31  had  been  obtained  for  a  restricted  Hartree-Fock  (HF)  reference  wave  function  with  a  distance-independent  spin-orbit  coupling  constant  according  to  Cohen 
and  Schneider  (denoted  as  “semiempirical”  in  the  table).  In  Ref.  31,  the  HF  orbitals  were  obtained  with  the  core  electrons  replaced  by  an  effective  core 
potential  (ECP)  and  including  a  core  polarization  potential. 


Method 

MRCI+Q(DKH2) 

CASPT2(DKH2) 

MRCI+Q(DKH2) 

HF-CCSD-T(ECP)‘ 

No.  of  correlated  e~ 

35 

35 

15 

35 

No.  of  configurations 

137  808  730 

6 

11  162  040 

basis  set 

{22s  1  \9p/  \  3d/5f/3g) 

[\0s/9p/Sd/5f/3g] 

{22s/\9pl  \3dl5fl3g) 

[6s/6p/6d/6fl2g] 

2K  (P"i) 

307.5 

303.4 

310.6 

o>,  (cnr1) 

136.1 

134.0 

124.3 

MRCI + Q(DKH2 + SO) 

CASPT2(DKH2  +  SO) 

MRCI + Q(DKH2 + SO) 

HF-CCSD-T(ECP+ SO)1 

Plus  SOC 

Breit-Pauli 

AMFI 

Breit-Pauli 

Semiempirical 

SOC  basis  set 

{22s/\9p/\3d) 

[\0s/9p/8d/5f/3g] 

(22sl\9pl\3d) 

I(l/2«):  Re  (pm) 

311.5 

306.5 

315.8 

311.4 

1(1  /2m):  a )e  (cm-1) 

124.5 

124.5 

112.6 

123.2 

“Reference  31. 


lence  region  (the  5s  and  5 p  shell)  of  Xe  (i.e.,  correlating  only 
15  electrons).  The  spin-orbit  coupled  electronic  states  were 
obtained  in  a  state-interaction  approach,  by  diagonalizing  the 
matrix  representation  of  the  electronic  Hamiltonian  and  the 
spin-orbit  Hamiltonian  in  a  basis  of  the  eigenfunctions  of  the 
electronic  Hamiltonian.  The  Breit-Pauli  Hamiltonian  was 
used  to  generate  the  spin-orbit  matrix  elements  for  the  inter¬ 
nal  configurations.  For  external  configurations  an  effective 
mean-field  Fock  operator  was  employed.62  The  spin-orbit 
matrix  was  constructed  in  a  smaller  (22s/  \9p/\3d)  basis  set, 
but  the  corresponding  diagonal  elements  were  shifted  in  en¬ 
ergy  to  match  the  Davidson-corrected  MRCI  energies.  An 
overview  of  selected  spectroscopic  parameters  of  the  A  22* 
and  I(l/2w)  states  and  the  approximations  involved  in  the 
quantum  chemical  calculations  is  given  in  Table  I. 

The  Davidson-corrected  MRCI  model  (MRCI+Q)  with 
35  correlated  electrons  features  the  best  agreement  with  the 
coupled-cluster  data  provided  by  Paidarova  and  Gadea.31 
The  equilibrium  bond  distance  and  the  harmonic  vibrational 
constant  are  remarkably  well  reproduced.  The  CASPT2  cal¬ 
culation  shows  larger  deviations  when  compared  to  Ref.  3 1 . 
The  equilibrium  bond  distance  Re  is  underestimated  by  about 
4  pm  in  the  CASPT2  +  SO  calculations.  It  is  overestimated 
by  the  same  amount  in  the  case  of  the  15-electron  MRCI 
+  Q+SO  results.  The  15-electron  MRCI +Q  wave  function  is 
not  a  good  approximation  to  the  35-electron  MRCI +  Q 
model  since  core  correlation  needs  to  be  taken  into  account 
explicitly. 

Concerning  the  DKH  order  we  note  that  the  equilibrium 
bond  lengths  of  the  22*  ground  state  changes  by  less  than 
0.5  pm  when  switching  from  the  DKH2  to  the  DKH8  Hamil¬ 
tonian  in  the  CASPT2  calculations.  Hence,  this  effect  can  be 
neglected  and  the  second-order  Hamiltonian  DKH2  is  suffi¬ 
ciently  reliable. 

III.  EXPERIMENT 

A  (1  +  1')  resonance-enhanced  two-photon  excitation 
scheme  was  used  to  excite  the  lowest  electronic  states  of  Xe£ 


via  selected  vibrational  levels  of  the  CO*  Rydberg  state  of 
Xe2  located  just  below  the  Xe!|,([5p]565'[l/2]1)  +  Xe(150) 
dissociation  limit.  The  vibrational  levels  of  the  C  0*  Rydberg 
state  were  chosen  as  intermediate  levels  because  they  are 
long  lived  and  have  been  well  characterized.63  64 

The  experimental  setup  used  in  the  present  study  has 
been  described  in  Refs.  65-67  and  consists  of  a  source  of 
coherent,  tunable  vacuum  ultraviolet  (VUV)  radiation  and  an 
electron/ion  time-of-flight  (TOF)  mass  spectrometer. 

The  laser  system  consists  of  three  dye  lasers  pumped  by 
the  doubled  and  tripled  output  of  a  Nd:YAG  (yttrium  alumi¬ 
num  garnet)  laser  at  a  repetition  rate  of  16^  Hz.  The  first  two 
dye  lasers  were  used  to  generate  tunable  VUV  radiation  by 
two-photon  resonance-enhanced  four-wave  mixing  in  kryp¬ 
ton.  The  output  of  the  first  dye  laser  was  tripled  using  two 
successive  /2-barium  borate  (BBO)  crystals  and  kept  fixed  so 
as  to  correspond  to  the  4p55p[l /2]0<-4p6(150)  two-photon 
resonance  in  krypton  [2^=94  092.86  cm"1  (Refs.  68  and 
69)].  The  fundamental  output  of  the  second  dye  laser  was 
overlapped  with  that  of  the  first  one  and  focused  into  a  gold- 
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FIG.  1.  PFI-ZEKE  photoelectron  spectrum  of  the  I(3/2g)  and  I(3/2w)  states 
of  ,32Xe,32Xe+  recorded  via  the  o' =21  vibrational  level  of  the  CO*  slate. 
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TABLE  II.  Measured  positions  Pobs  and  differences  between  measured  and  calculated  positions  (Ay=yobs 
-i'cak)  of  the  vibrational  levels  of  the  I(3/2g)  state  of  l3lXe,32Xe+  relative  to  the  X  0*  (u"=0)  ground  neutral 
state.  The  uncertainties  include  possible  errors  in  the  calibration  and  in  the  determination  of  the  field-induced 
shift  of  the  ionization  thresholds.  DL  indicates  the  position  of  the  Xe+(2/>3/2)+Xe(l50)  dissociation  limit. 


vobs  (cm  *) 

Ay  (cm  *)* 

v+ 

vobs  (cm-1) 

AP  (cm-1)* 

v+ 

vobi  (cm  ’) 

Ay  (cm-1)* 

0 

96  219.9(20)b 

0.3 

18 

97  102.5 

36 

97  658.7(17) 

-0.02 

1 

96  278.0(3  l)b 

0.6 

19 

97  141.9 

37 

97  681.2(18) 

0.7 

2 

96  334.7(36)b 

0.6 

20 

97  179.3(ll)b 

-0.9 

38 

97  701.3(25) 

-0.05 

3 

96  389.9(27)b 

0.05 

21 

97  217.3(38 )b 

-0.2 

39 

97  721.3(17) 

-0.01 

4 

96  445.0(22)b 

0.4 

22 

97  253.8 

40 

97  740.0(20) 

-0.4 

5 

96  497.7(39)b 

-0.5 

23 

97  289.1 

41 

97  758.2(16) 

-0.4 

6 

96  550.4(47)b 

-0.5 

24 

97  323.4 

42 

97  775.4(16) 

-0.5 

7 

96  602.2(26)b 

-0.3 

25 

97  357.2(19) 

0.5 

43 

97  792.0(16) 

-0.4 

8 

%  652.4(38)b 

-0.7 

26 

97  389.3(19) 

0.2 

44 

97  807.2(17) 

-0.9 

9 

96  701.7(1  l)b 

-0.9 

27 

97  420.9(15) 

0.5 

45 

97  822.6(16) 

-0.4 

10 

96  750. 1(l8)b 

-1.0 

28 

97  450.9(17) 

0.2 

46 

97  836.3(29) 

-0.8 

11 

96  797.0(l7)b 

-1.6 

29 

97  480.3(14) 

0.2 

47 

97  849.3(25) 

-1.1 

12 

96  844.2(25)b 

-0.9 

30 

97  508.8(15) 

0.3 

48 

97  861.5(28) 

-1.4 

13 

96  889.4(  1 1  )b 

-1.2 

31 

97  535.9(24) 

0.01 

49 

97  873.7(17) 

-LI 

14 

96  935.0(l0)b 

0.01 

32 

97  562.5(22) 

0.1 

50 

97  885.1(20) 

-0.8 

15 

96  978.9(30)b 

0.5 

33 

97  588.4(14) 

0.5 

51 

97  895.9(10) 

-0.5 

16 

97  021.2(28)b 

0.4 

34 

97  612.9(18) 

0.5 

52 

97  905.1(16) 

-LI 

17 

97  062.3(10)b 

0.1 

35 

97  636.5(17) 

0.5 

DL 

98  019.78 

‘‘When  the  transitions  could  not  be  observed,  the  calculated  positions  are  listed, 
^ata  from  Rupper  et  al.  (Ref.  32). 


coated  four-wave  mixing  cell  with  a  slow  constant  flow  of 
krypton  at  a  pressure  of  —16  mbars  to  generate  the  VUV 
radiation  required  to  drive  the  CO*  <—  X0+  transition  of  Xe2. 
The  VUV  radiation  was  separated  from  the  fundamental  la¬ 
ser  beams  by  a  MgF2  prism  and  directed  toward  the  photo¬ 
excitation  region  where  it  intersected  the  probe  gas  beam  and 
induced  a  transition  to  a  selected  vibrational  level  of  the  C  0* 
state.  The  output  of  the  third  dye  laser  was  used  to  access  the 
electronic  states  of  Xe2  from  the  selected  intermediate  levels. 
The  calibration  of  the  fundamental  wave  number  of  the  third 
dye  laser  was  achieved  by  splitting  off  a  small  fraction  of  the 
laser  beam  and  directing  it  toward  an  optogalvanic  (OG)  cell 
filled  with  neon  or  argon  and,  simultaneously,  recording  an 
etalon  spectrum. 

Xe2  rare  gas  dimers  were  generated  in  a  pulsed  super¬ 
sonic  expansion  of  pure  xenon  (stagnation  pressure  of 
—3.6  bars)  through  a  pulsed  solenoid  valve  (orifice  diameter 
of  0.4  mm).  Before  entering  the  photoexcitation  region  the 


supersonic  beam  passed  through  a  skimmer  (orifice  diameter 
of  1.0  mm).  The  supersonic  beam  crossed  the  laser  beams  at 
right  angles. 

PFI-ZEKE  photoelectron  spectra  of  Xe2  were  recorded 
by  detecting  the  electrons  produced  by  delayed  pulsed  field 
ionization  of  very  high  Rydberg  states  located  below  the  vi- 
bronic  states  of  Xe2  using  a  sequence  of  two  successive  elec¬ 
tric  field  pulses  as  a  function  of  the  wave  number  F3  of  the 
fundamental  or  doubled  output  of  the  third  dye  laser. 

IV.  RESULTS 

A.  The  \(3/2g)  and  l(3/2u)  states  of  Xe£ 

PFI-ZEKE  photoelectron  spectra  following  single¬ 
photon  excitation  from  the  neutral  ground  state  provided  in¬ 
formation  only  on  the  low  vibrational  levels  of  the  1(3  !2g) 
and  I(3/2m)  states  of  Xe2  up  to  u+=21  and  u+=2, 
respectively.32  Close  to  the  Xe(150)  +  Xe+(2P3/2)  dissociation 


TABLE  III.  Measured  positions  yobs  and  differences  between  measured  and  calculated  positions  (Ay=yobs 
-ycajc)  °f  the  vibrational  levels  of  the  I(3/2m)  state  of  ,3lXe,32Xe+  relative  to  the  X0*  (i/'=0)  ground  neutral 
state.  The  uncertainties  include  possible  errors  in  the  calibration  and  in  the  determination  of  the  field-induced 
shift  of  the  ionization  thresholds.  DL  indicates  the  position  of  the  Xe+(2/>3/2)+Xe(,50)  dissociation  limit. 


S'obs  (cm  ‘) 

Ay  (cm  ,)a 

u+ 

yobs  (cm-1) 

Ay  (cm"1)1 

u+ 

yob8  (cm"1) 

Ay  (cm"1)1 

0 

97  577.6(15) 

0.3 

8 

97  734.1 

16 

97  846.8(17) 

-0.1 

1 

97  600.0(15) 

0.5 

9 

97  750.5 

17 

97  858.1(16) 

-0.2 

2 

97  621.4(16) 

0.4 

10 

97  766.5(19) 

0.3 

18 

97  868.8(15) 

-0.3 

3 

97  642.2(12) 

0.5 

11 

97  781.6(17) 

0.4 

19 

97  879.0(24) 

-0.3 

4 

97  662.2(16) 

0.6 

12 

97  795.7(15) 

0.1 

20 

97  888.8(16) 

-0.3 

5 

97  681.2(19) 

0.4 

13 

97  809.5(28) 

0.2 

21 

97  896.9(13) 

-1.4 

6 

97  699.8(23) 

0.5 

14 

97  822.6(16) 

0.2 

22 

97  907.0(12) 

-0.1 

7 

97  717.6(15) 

0.5 

15 

97  835.2(27) 

0.2 

DL 

98  019.78 

“When  the  transitions  could  not  be  observed,  the  calculated  positions  are  listed. 
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FIG.  2.  PFI-ZEKE  photoelectron  spectra  of  the  1(3/2 g)  and  1(3 /2m)  states 
just  below  the  Xe+(2P3/2)  +  Xe(,50)  dissociation  limit  of  ,29Xe,32Xe+  re¬ 
corded  via  the  u'=20  vibrational  level  of  the  CO*  state  (upper  trace)  and  of 
,3,Xe,32Xe+  recorded  via  the  i/=21  and  22  vibrational  levels  of  the  CO* 
state  (lower  traces). 


limit  the  interpretation  of  the  single-photon  PFI-ZEKE  pho¬ 
toelectron  spectrum  was  complicated  by  atomic  lines  belong¬ 
ing  to  the  Rydberg  series  converging  to  the  2P3/2  dissociation 
limit  of  atomic  xenon.32  These  undesirable  atomic  contribu¬ 
tions  could  be  avoided  in  the  present  work  using  a 
resonant — and  thus  species-selective — two-photon  excitation 
sequence.  Moreover,  different  regions  of  the  ionic  potential 
energy  functions  could  be  accessed  by  selecting  different  vi¬ 
brational  levels  of  the  intermediate  state. 

The  spectrum  of  the  1(3/2 g)  and  1(3  /2m)  states  of 
131Xe132Xe*  recorded  via  the  v'  =  21  vibrational  level  of  the 
CO*  state  ls  presented  in  Fig.  1.  Long  vibrational  progres¬ 
sions  were  observed  and  assigned  to  levels  as  high  as  i/* 
=  52  for  the  1(3/2 g)  state  and  t/*=22  for  the  1(3  /2m)  state. 
The  unambiguous  vibrational  assignments  were  derived  from 
an  analysis  of  the  isotopic  shifts  using  the  procedure  de¬ 
scribed  in  Refs.  70  and  71.  In  the  wave  number  region  be¬ 
yond  97  910  cm-1,  an  unambiguous  assignment  of  the  spec¬ 
tral  features  could  not  be  derived  because  of  (1)  the  spectral 
congestion,  (2)  the  weakness  of  most  spectral  features  and 
(3)  the  irregularity  of  the  intensity  distribution  partially 
caused  by  channel  interactions  with  autoionizing  Rydberg 
states.72  These  autoionizing  states  were  independently  ob¬ 
served  in  photoionization  spectra  recorded  via  the  same  in¬ 
termediate  vibrational  levels  of  the  CO*  state  (not  shown 
here).  Spectral  lines  in  the  photoelectron  spectra  observed  in 
this  region  could  correspond  to  vibrational  levels  of  any  of 
the  four  ionic  states  of  Xe2  associated  with  the  Xe+(2P3/2) 
+  Xe(1S0)  dissociation  limit.  The  positions  of  the  vibrational 
levels  for  the  I(3/2g)  and  I(3/2m)  states  of  131Xe132Xe*  are 
summarized  in  Tables  II  and  III. 

The  PFI-ZEKE  photoelectron  spectra  of  the  1(3/2 g)  and 
1(3/ 2m)  states  of  129Xe132Xe*  and  131Xe132Xe*  between 
97  820  and  98  000  cm"1  are  compared  in  Fig.  2.  The  inten¬ 
sity  distributions  strongly  depend  on  the  selected  vibrational 
level  of  the  intermediate  CO*  state  because  of  different 
Franck-Condon  factors  and  channel  interactions.72  The  top 
trace  in  Fig.  2  corresponds  to  the  spectrum  of  the  129Xe132Xe 


isotopomer  recorded  via  the  C  0*  (i/  =  20)  intermediate  level 
whereas  the  bottom  two  traces  correspond  to  spectra  of 
131Xe132Xe  recorded  through  the  C  0*  (i/  =  21  and  22)  levels, 
as  indicated.  The  line  marked  by  an  asterisk  in  the  lower  two 
traces  of  Fig.  2  could  not  be  assigned  to  either  the  1(3/2 g)  or 
the  1(3/2 u)  states  of  Xe2.  This  line  could  correspond  to  a  low 
vibrational  level  of  the  1(1/2 g)  state  or  to  a  high  vibrational 
level  (i/*>110)  of  the  1(1 /2m)  state. 

B.  The  11(1 /2u)  state  of  XeJ 

In  our  previous  study  of  the  lowest  electronic  states  of 
Xe2  following  single-photon  excitation  from  the  X  0*  ground 
state  of  Xe2,  only  transitions  to  the  lowest  twelve  vibrational 
levels  of  the  11(1 /2m)  state  could  be  observed.32  An  unam¬ 
biguous  assignment  of  the  absolute  value  of  the  vibrational 
quantum  number  was,  however,  hindered  by  the  fact  that  the 
isotopic  shifts  were  not  large  enough. 

The  isotope  selectivity  of  the  resonant  two-photon  exci¬ 
tation  sequence  used  in  the  present  investigation  allowed  the 
observation  of  the  11(1  /2m)(u+)<—  X 0*  progression  for  dif¬ 
ferent  isotopomers  of  Xe2.  The  PFI-ZEKE  photoelectron 
spectrum  of  the  11(1 /2m)  state  of  131Xe132Xe*  recorded  via 
the  l?' =21  vibrational  level  of  the  CO*  state  of  Xe2  is  dis¬ 
played  in  the  lower  panel  of  Fig.  3  and  the  corresponding 
line  positions  are  listed  in  Table  IV.  The  decreasing  back¬ 
ground  signal  in  the  photoelectron  signal  between  107  100 
and  107  350  cm"1  was  caused  by  a  strong  electron  signal 
close  to  the  time  gate  of  the  background  window  in  the  TOF 
trace.  These  prompt  electrons  are  mostly  generated  by  pro¬ 
cesses  where  the  laser  beams  hit  metallic  surfaces.  The  line 
marked  by  an  asterisk  corresponds  to  an  impurity  line.  From 
the  analysis  of  the  isotopic  shifts  of  the  sixth  line  in  the 
vibrational  progression  [see  upper  panels  in  Fig.  3  for  a  com¬ 
parison  of  the  assignments  to  i/*=5  (full  lines)  and  to  i/*=6 
(dashed  lines)],  the  vibrational  assignment  of  the  11(1 /2m) 
state  tentatively  proposed  by  Rupper  et  m/.32  could  be 
confirmed. 

Unfortunately,  the  11(1  Hg)  state  could  not  be  observed 
from  any  of  the  selected  vibrational  levels  of  the  intermedi¬ 
ate  C  0*  state  and  we  are  not  able  to  confirm  or  refute  the 
assignment  of  the  single  line  observed  in  our  previous 
study.32 

C.  Spectroscopic  constants  of  the  l(3/2g)  and  l(3/2u) 
states  of  Xe2 

Analyzing  the  vibrational  progressions  associated  with 
the  I(3/2g)  and  I(3/2m)  states  of  131Xe132Xe2  in  terms  of  the 
standard  expansion  formula 

v  =  EJhc  +  cu*(i/*  +  1/2)  -  ovt*(i/*  +  1/2)2 

+  w,y>+  +  1/2)3  -  K/2  -  (OfX+14  +  «,£/8]  ( 1 ) 

led  to  the  set  of  adiabatic  ionization  energies  harmonic 
(g>*)  and  anharmonic  (ov**,  <*>*>£,  cu^z*)  vibrational  constants 
summarized  in  Table  V.  Spectroscopic  constants  for  the 
1(3/2 g)  state  were  derived  using  the  transition  wave  numbers 
determined  from  the  spectra  shown  in  Figs.  1  and  2  and  the 
positions  of  the  lower  vibrational  levels  reported  by  Rupper 
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FIG.  3.  Upper  panels:  upper  traces:  PFI-ZEKE  photo¬ 
electron  spectra  of  the  11(1  Hu)  (u+=5)  state  of  differ¬ 
ent  isotopomers  of  Xe2  recorded  via  the  u'=21  vibra¬ 
tional  level  of  the  CO*  state.  Lower  inverted  traces: 
Simulated  spectra  assuming  that  the  ionic  state  has  vi¬ 
brational  quantum  number  u+=5  (full  line)  and  u+=6 
(dashed  line):  (a)  l29Xel29Xe,  (b)  l29Xe132Xe,  (c) 
l32Xe132Xe,  and  (d)  ,32Xe136Xe.  Lower  panel:  PFI- 
ZEKE  photoelectron  spectrum  of  the  11(1  Hu)  state  of 
,32Xe,32Xe*  recorded  via  the  v'  =  21  vibrational  level  of 
the  CO*  state.  The  line  marked  by  an  asterisk  corre¬ 
sponds  to  an  impurity  line. 


et  al 32  The  constants  of  the  1(3  /  2w)  state  were  derived  solely 
from  the  experimental  data  recorded  in  this  work.  In  addi¬ 
tion,  values  for  the  dissociation  energies  Dq  were  determined 
from 

D+0(Xep  =  £,(Xe,2P,+  -  lS0)lhc 

-  EiiXcj  <-  Xe2,X  0p/hc  +  D0(Xe2,X  Op  (2) 

using  the  atomic  ionization  energies  of  97  833.783  cm-1  cor¬ 
responding  to  the  formation  of  the  2P3/2  state  of  Xe+  (Ref. 
73)  and  the  dissociation  energy  D0=186cm_1  of  the  X0* 
ground  state  of  Xe2M  The  spectroscopic  constants  summa¬ 
rized  in  Table  V  are  consistent  with  the  earlier  results  of 
Tonkyn  et  al.,39  Lu  et  al.3*  and  Rupper  et  al  31  but  represent 
a  more  complete  set. 

V.  THE  POTENTIAL  ENERGY  FUNCTIONS  OF  THE  SIX 
LOW-LYING  ELECTRONIC  STATES  OF  Xe+ 

A.  Potential  model 

The  semiempirical  potential  model  first  presented  in  Ref. 
29  for  the  determination  of  the  six  lowest  electronic  potential 
functions  of  Ar2  and  later  used  for  Kr2  (Ref.  30)  and  Xe2 


(Ref.  32)  was  based,  as  all  previous  ab  initio  studies,31,40-43 
on  the  approximation  of  an  /^-independent  spin-orbit  cou¬ 
pling  constant  introduced  by  Cohen  and  Schneider.44  The 
consequences  of  this  approximation  have  never  been  quanti¬ 
fied  so  far.  Relativistic  ab  initio  quantum  chemical  calcula¬ 
tions  (see  Sec.  II)  allowed  the  investigation  of  the 
R-dependence  of  the  spin-orbit  coupling  constant  a(R)  by 
comparing  the  potential  energy  functions  of  the  Ilw  and  n^, 
states  derived  in  scalar-relativistic  calculations  with  those  for 
the  1(3/2 u)  and  1(3/2 g)  states  determined  in  relativistic  cal¬ 
culations  including  also  spin-orbit  coupling.  According  to  the 
spin-orbit  interaction  matrix  presented  in  Table  VI  the  poten¬ 
tial  energy  functions  of  the  2II  state  and  the  2ll3/2  state  differ 
by  a(R)/2.  Thus  a(R)  can  be  determined  at  each  intemuclear 
separation  for  which  both  calculations  have  been  performed 
by  simple  subtraction. 

The  upper  sets  of  open  squares  and  circles  in  Fig.  4 
correspond  to  the  values  of  a(R)  determined  directly  from 
the  MRCI+2  (DKH2  +  SO)  ab  initio  calculations  (see  col¬ 
umn  1  in  Table  I).  At  large  distances,  the  calculated  spin- 
orbit  coupling  constants  tend  asymptotically  to  a  value  larger 
than  the  atomic  value,  both  for  the  g  and  the  u  states.  To 


TABLE  IV.  Measured  positions  vobs  and  differences  between  measured  and  calculated  positions  (AF=Fobs 
-Fcak;)  of  the  vibrational  levels  of  the  11(1  /2m)  and  11(1 12g)  states  of  ,3,Xe,32Xe*  relative  to  the  X0*  (i/'=0) 
ground  neutral  state.  The  uncertainties  include  possible  errors  in  the  calibration  and  in  the  determination  of  the 
field-induced  shift  of  the  ionization  thresholds.  DL  indicates  the  position  of  the  Xe+(2/>1/2)  +  Xe(,50)  dissociation 
limit. 

y+ 

Fobs  (cm"1) 

A?  (cm-1) 

v * 

Fobs  (cm-1) 

AF  (cm  ’) 

0 

107  159.6(23) 

-1.6 

8 

107  517.8(29) 

-1.2 

1 

107  207.7(17) 

-1.4 

9 

107  558.7(38) 

-1.0 

2 

107  254.5(24) 

-1.6 

10 

107  603.2(29)a 

3.8 

3 

107  301.2(19) 

-1.0 

11 

107  643.9(16)“ 

5.7 

4 

107  346.4(17) 

-1.0 

5 

107  391.0(19) 

-0.7 

11(1 /2g) 

6 

107  434.3(19) 

-0.7 

0 

108  132.3(13)a,b 

7 

107  476.8(15) 

-0.7 

DL 

108  556.71 

“Data  from  Rupper  et  al.  (Ref.  32),  incorrectly  assigned  to  the  l3,Xel32Xe+  isotopomer.  The  corrected  positions 

and  deviations  are  107  598.5(10)  and  -0.9  cm-1  for  u+=  10  and  107  637.0(5)  and  -1.2  cm-1  for  v+=  1 1,  respec¬ 
tively. 

‘"Data  from  Rupper  et  al.  (Ref.  32).  The  values  corresponds  to  131 3Xe2. 
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234306-7  fl-dependence  of  spin-orbit  coupling  in  XeJ 


TABLE  V.  Adiabatic  ionization  energies  Eh  dissociation  energies  Dj,  and  vibrational  constants  a>*,  (oey+ 
of  the  I(3/2g)  and  I(3/2u)  states  of  ,3,Xe,32Xe+. 


State 

£,  {he  cm-1) 

Dq  (cm-1) 

o>*  (cm*1) 

<*>Xt  (cm-1) 

cD,y*  (cm"1) 

Reference 

1(3/2*) 

%  220.7(15) 

1799.1(15) 

58.17(35) 

0.456(26) 

-0.001  98(73)* 

This  workb_d 

96  220.2(10) 

1799.6(15) 

58.61(27) 

0.506(13) 

32b-d 

96  226(4) 

1793(5) 

58.36(34) 

0.484(15) 

38" 

%  359(24) 

1670(24) 

34" 

>1220(5)' 

>45' 

39" 

96  327(40) 

1694(40) 

35" 

%  545(121) 

1492(121) 

33d-* 

96  343(40) 

1815(40) 

55.5(28) 

81" 

1565(40) 

55(4) 

82" 

%  143 

1876 

58.74 

0.508 

3l" 

1178 

49.7 

0.3 

40" 

I(3/2«) 

97  577.6(15) 

442.2(15) 

23.19(54) 

0.419(54) 

0.002  7(15) 

This  workb_d 

97  576.6(10) 

443.2(15) 

21.89(40) 

32b_c 

97  582(4) 

437(5) 

38" 

97  214(24) 

807(24) 

34" 

97  576(2) 

442(2) 

23.1 

0.55 

39" 

97  617(40) 

403(40) 

35" 

<403 

33" 

97  512(40) 

645(40) 

33.0(17) 

81" 

613(40) 

33(4) 

82" 

97  437 

582 

26.53 

0.425 

31" 

306 

58.5 

0.3 

40" 

aw^;=-0.000  026  6(68)  cnT1. 

^The  uncertainties  in  £,  and  DjJ  include  the  full  width  at  half  maximum  of  the  observed  transitions  and  potential 
errors  in  the  wave  number  calibration  and  in  the  determination  of  the  field-induced  shifts  of  the  ionization 
thresholds. 

cThe  uncertainty  in  the  vibrational  constants  represent  one  standard  deviation  in  the  fit. 
dExperimental  values. 

eValues  correspond  to  a  mixture  of  all  isotopomers  which  is  best  represented  by  a  Active  average  isotopomer  of 
Xe2  with  reduced  mass  /u= 65.646  u. 

fTonkyn  and  White  (Ref.  39)  did  not  derive  an  absolute  vibrational  assignment  of  the  1(3 /2g)  state. 
gSemiempirical  values  using  known  spectroscopic  data. 

‘'Theoretical  values. 


correct  for  this  artefact  of  the  calculations,  the  curves  were 
scaled  by  a  factor  of  0.933  for  the  u  states  and  0.928  for  the 
g  states  so  that  the  asymptotic  value  of  a(R)  exactly  corre¬ 
sponds  to  the  spin-orbit  splitting  of  the  Xe+  2P  ground  state. 
The  scaled  ag/u(R)  values  are  also  displayed  in  Fig.  4  (lower 
sets  of  open  squares  and  circles).  The  spin-orbit  coupling 
constant  varies  by  about  10%  over  the  range  of  intemuclear 
distance  relevant  for  the  present  study.  Consequently,  the 
^-dependence  of  the  spin-orbit  coupling  constant  must  be 
included  in  a  quantitative  description  of  the  potential  energy 
functions  of  the  six  low-lying  electronic  states  of  XeJ.  We 


TABLE  VI.  Spin-orbit  interaction  matrix  in  Hund’s  case  (a)  basis  describing 
the  coupling  between  the  states  of  2£*  and  2I1  symmetry  in  the  homonuclear 
rare  gas  dimer  ions. 


2T+ 

^  1/2 

2n„2 

2n  „ 

2s;,2 

ViW 

a(R) 

~  Vi 

2n,« 

a(R) 

~  Vi 

a(R) 

VUW+  — 

a(R) 

V„(R)~— 

should  note  that  also  the  minimal  CASPT2  model  (see  col¬ 
umn  2  in  Table  I)  features  a  very  similiar  behavior  of  a(R) 
although  these  a(R)  curves  are  below  the  asymptotic  limit 
and  would  require  a  scaling  factor  that  is  slightly  larger  than 
1 .  However,  in  view  of  the  validation  of  the  MRCI  data  in 
Sec.  II,  we  consider  only  the  MRCI +2  (DKH2  +  SO)  data  in 
the  following. 

We  found  that  the  R-dependence  of  the  scaled  spin-orbit 
coupling  constant  a{R)  of  the  u  and  g  states  could  be  de¬ 
scribed  almost  perfectly  by  two  rather  simple  analytical  ex¬ 
pressions, 

au(/e)  =  aatom-au0(l  -(1  (3) 

ag(R)  =  aMom  +  ag0e-atiR°,\  (4) 

where 

tfa«on>  =  aXe  =  jAXe  =  7024.61 7  cm'1 

(^xe=[^(Xe+,2P1/2)-£(Xe*,2/>j/2)]//»c).74  Equation  (3)  is  a 
Morse-type  function  and  Eq.  (4)  falls  exponentially  at  in¬ 
creasing  R.  The  optimal  parameters  aui  and  a#l(i=0, 1 ,2) 
are  listed  in  the  lower  part  of  Table  VII  and  the  correspond- 
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FIG.  4.  /^-dependent  spin-orbit  coupling  constant  a{R)  for  the  u  (circles) 
and  the  g  (squares)  states  derived  from  the  comparison  of  the  relativistic 
(including  spin-orbit  coupling)  with  the  scalar-relativistic  potential  energy 
curves  derived  by  MRCI+Q  (DKH2)  ab  initio  calculations  (see  column  1  in 
Table  I)  (upper  pair  of  curves).  The  lower  pair  of  curves  are  scaled  such  that 
the  asymptotic  value  exactly  corresponds  to  the  spin-orbit  splitting  of  the 
Xe+  2P  ground  state.  The  full  lines  represent  fits  based  on  the  expressions  in 
Eqs.  (3)  and  (4). 

ing  au{R)  and  ag{R)  curves  are  displayed  as  full  lines  in 
Fig.  4. 

The  effects  of  the  spin-orbit  interaction  are  treated  using 
the  same  coupling  matrix  as  proposed  by  Cohen  and 
Schneider,44  improved  by  the  /^-dependence  of  the  spin-orbit 
coupling  constant.  The  procedure  consists  of  diagonalizing, 
at  each  value  of  the  intemuclear  distance  R,  the  spin-orbit 
interaction  matrix  given  in  Table  VI.  Because  the  spin-orbit 
operator  only  couples  states  of  the  same  g/u  symmetry,  two 
distinct  3X3  matrices  must  be  set  up,  one  for  the  g  states, 
the  other  for  the  u  states.  The  spin-orbit  interaction  matrix  is 
expressed  in  the  basis  set  adequate  for  the  short-range  part  of 
the  potentials  in  which  the  quantum  numbers  A  and  2  cor¬ 
responding  to  the  projection  of  the  total  electronic  orbital  and 
spin  angular  momenta  onto  the  intemuclear  axis  are  good 
quantum  numbers.  Consequently,  the  electronic  potential  en¬ 


ergies  of  the  2£|/2, 2H1/2,  and  2n3/2  states  appear  as  diagonal 
elements.  The  spin-orbit  operator  is  not  only  described  by 
diagonal  contributions  of  ±a(R)/2  which  induce  a  splitting 
of  the  2n  state  into  the  two  components  2ni/2  and  2n3/2,  but 
also  by  off-diagonal  elements  -a{R)/\^2  which  couple  states 
of  the  same  value  of  O  (the  projection  of  the  total  angular 
momentum  onto  the  intemuclear  axis),  namely,  the  2X|/2  and 
2ni/2  states. 

The  potential  energy  functions  of  the  22+  and  2 II  states 
under  neglect  of  the  spin-orbit  interaction  are  expressed  as 

V^R)  =  Axe~b^R  - 

3  c 

-2 /*,(*,*  aHSt  +  V*-  (5) 

n= 2  K 

with  A  =  2,n.  Because  the  22*  state  is  only  weakly  bound, 
the  second  term  in  Eq.  (5),  which  describes  the  chemical 
bond,  is  neglected  for  this  state  ( B %  =0).  Vdiss  is  a  constant 
used  to  relate  the  potential  energies  to  the  energy  of  the  X  0* 
(i/'=0)  ground  neutral  state  and  is  defined  as 

Vdiss  =  Oo(Xe2,X  Op  +  E,(Xe,2PV2 «-  lS0)/hc  + 

(6) 

with  aatom=aXc.  The  value  of  Vdiss=101  532.1  cm-1  is  ob¬ 
tained  from  the  literature  values  ^atom=3^xe 
=7024.617  cm'1,74  E^Xe^P^)/ hc=91  833.783  cm"1/5 

and  D0(Xe2,X 0*)=  186  cm’1.64 

Because  the  and  2n  potential  energy  curves  are  cor¬ 
related  with  the  Xe(150)-»-Xe+(2P)  dissociation  limit,  the 
long-range  behavior  of  these  states  is  described  in  terms  of 
the  interaction  between  a  neutral  Xe  atom  and  a  Xe+  ion.  The 
long-range  coefficients  are  equal  for  the  u  and  the  g  states. 
Values  for  the  C4  and  C6  coefficients  of  Xe2  were  derived 
using  the  formalism  described  in  the  appendix  of  Ref.  67. 

The  potential  energy  functions  Vj(R)  of  the  three  u  and 
the  three  g  states  were  then  calculated  from  the  Vj  (/?), 


TABLE  VII.  Parameters  of  the  interaction  potentials  of  the  lowest  electronic  states  of  Xe^.  The  uncertainties 
represent  95%  confidence  intervals.  Parameters  given  without  uncertainties  were  held  fixed. 


u  states 

g  states 

2s; 

2n„ 

2X+ 

K 

2ns 

*.,A  (A) 

3.077  89  ±0.000  30 

4.607  05  ±  0.000  51 

7.5* 

3.6960  ±  0.0022 

D.,a  (cm-1) 

10  634.87  ±0.62 

470.967  ±  0.088 

1772.25  ±0.82 

P\ 

1.290  29  ±0.000  47 

4.0204  ±0.0024 

1.5162  ±  0.0024 

frA  (A-1) 

2.151  50±  0.00041 

2.842  95  ±0.000  53 

1 . 1 1 57  ±  0.00 1 0  2.6963  ±0.0019 

c4.  A  (cm-1  ATb 

234  837 

234  837 

234  837 

234  837 

c6.  A  (cm-1  A6)*-” 

2  331  794 

1  167  755 

2  331  794 

1  167  755 

A\  (cm-1)' 

2.458  X  107 

1.094  X  108 

148  617 

4.858  X  107 

B  \  (cm-1)' 

6.764  X  106 

1471 

1.743  XI 06 

au(R) 

a0  (cm-1) 

597.2* 

1084.5“ 

a,  (A-1) 

1.566* 

0.086“ 

«2 

2.359  A’ 

2.500“ 

“This  parameter  was  kept  constant  during  the  fit. 

bUsing  1£A=219  474.6  cm"1,  la0=0.529  177  21 1  A.  See  Appendix  in  Ref.  67. 

cDetermined  from  the  conditions  [dV\(J?)/dJ?]|*=i? A  =  0  and  V\(Re  A)  =  Vdlss-De  A.  See  Appendix  A  in  Ref.  29. 
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TABLE  VIII.  Measured  positions  Mobs  from  Rupper  et  al.  (Ref.  32)  and  differences  between  measured  and 
calculated  positions  (AM=Mobs- Mcak)  of  the  vibrational  levels  of  the  1(1 /2m)  state  of  ,3,Xe,32Xe*  relative  to  the 
X0+  (u"=0)  ground  neutral  state.  The  uncertainties  include  possible  errors  in  the  calibration  and  in  the  deter¬ 
mination  of  the  field-induced  shift  of  the  ionization  thresholds.  DL  indicates  the  position  of  the  Xe+^P^) 
+  Xe(,50)  dissociation  limit. 


V + 

Mobs  (cm-1)  Am  (cm"1)' 

v * 

Mo*  (cm-1)  Am  (era-1)1 

Mobs  (cm-1) 

Am  (cm*1)1 

0 

90  147.5(17) 

3.5 

21 

92  512.3(30) 

-2.3 

42 

94  455.2(12) 

-1.2 

I 

90  269.8(20) 

3.3 

22 

92  614.6(25) 

-2.2 

43 

94  537.5(23) 

-0.6 

2 

90  390.6(35) 

2.6 

23 

92  715.8(30) 

-2.3 

44 

94  618.4(15) 

-0.4 

3 

90  510.5(17) 

1.9 

24 

92  816.3(21) 

-2.1 

45 

94  697.9(13) 

-0.6 

4 

90  629.9(30) 

1.6 

25 

92  915.4(14) 

-2.3 

46 

94  777.2(23) 

-0.06 

5 

90  748.1(14) 

1.2 

26 

93  013.6(23) 

-2.5 

47 

94  855.5(17) 

0.5 

6 

90  866.0(13) 

1.4 

27 

93  110.7(22) 

-2.7 

48 

94  932.2(18) 

0.4 

7 

90  982.5(11) 

1.1 

28 

93  207.7(28) 

-2.1 

49 

95  008.4(15) 

0.8 

8 

91  097.4(19) 

0.2 

29 

93  302.8(30) 

-2.4 

50 

95  083.7(20) 

1.3 

9 

91  212.0 

30 

93  396.9(13) 

-2.7 

51 

95  157.2(19) 

1.0 

10 

91  325.2(30) 

-0.7 

31 

93  490.5(23) 

-2.6 

52 

95  230.6(15) 

1.5 

11 

91  438.1(24) 

-0.7 

32 

93  583.3(17) 

-2.3 

53 

95  303.3(16) 

2.4 

12 

91  550.2(19) 

-0.5 

33 

93  674.9(20) 

-2.2 

54 

95  374.2(18) 

2.4 

13 

91  660.6(28) 

-1.1 

34 

93  767.6 

55 

95  443.9(18) 

2.2 

14 

91  771.2(28) 

-0.5 

35 

93  855.1(14) 

-2.0 

56 

95  513.9(31) 

3.3 

15 

91  879.5(22) 

-1.2 

36 

93  944.0(14) 

-1.7 

57 

95  581.7(20) 

3.1 

16 

91  987.5(28) 

-1.3 

37 

94  031.8(33) 

-1.5 

58 

95  649.6(23) 

4.1 

17 

92  094.1(28) 

-1.8 

38 

94  118.2(26) 

-1.7 

59 

95  715.1(17) 

3.6 

18 

92  202.0 

39 

94  204.3(14) 

-1.2 

60 

95  781.0(14) 

4.5 

19 

92  305.7(30) 

-1.5 

40 

94  288.7(14) 

-1.4 

61 

95  845.5(37) 

5.0 

20 

92  409.8(29) 

-1.6 

41 

94  372.5(17) 

-1.3 

DL 

98  019.78 

“When  the  transitions  could  not  be  observed,  the  calculated  positions  are  listed. 


Vn  (/?),  (R),  and  Vn  (R)  potentials  in  a  straightforward 

manner  using  the  interaction  matrix  presented  in  Table  VI. 
The  calculated  transition  wave  numbers  were  determined  by 
numerically  solving  the  radial  Schrodinger  equation 


h2  d2 

■ - r  +  hcVXR) 

2m,  d/?2  J 


i/iiJR)  =  Eiju>fiiju(R) 


(7) 


corresponding  to  the  vibrational  motion  of  the  isotopomer 
with  reduced  mass  m,  in  the  potential  V}(/?)  (see  below).  In 
Eq.  (7),  the  indices  i,  j,  and  v  designate  the  isotopomer,  the 
electronic  state,  and  the  vibrational  quantum  number  of  the 
ion,  respectively. 


B.  Fitting  procedure 

The  parameters  of  the  model  potential  energy  functions 
[Eq.  (5)]  were  derived  from  the  experimental  data  in  a  non¬ 
linear  least-squares  fitting  procedure.  The  numerical  proce¬ 
dure  relied  on  a  discrete  variable  representation  (based  on  the 
Gauss-Chebyshev  quadrature  of  the  first  kind75-78)  on  a  grid 
with  801  equidistant  grid  points  in  the  range  of  intemuclear 
distances  /?=2.0- 10.0  A.  In  the  fitting  procedure,  all  param¬ 
eters  describing  the  long-range  interaction  were  kept  fixed  at 
the  values  listed  in  Table  VII. 


1.  u  states 

The  u  states  are  experimentally  well  characterized  and 
the  observed  and  unambiguously  assigned  vibrational  bands 
[59  bands  for  the  1(1  /2m)  state,  21  bands  for  the  1(3  /2m) 
state,  and  12  bands  for  the  11(1  /2m)  state]  form  a  large  data 
set  with  which  the  potential  energy  functions  of  the  three 


electronic  states  of  u  symmetry  can  be  determined.  A  total  of 
eight  parameters  (ReX,  De A,  ($ A,  fcA;  A=2,II)  were  refined 
in  the  fitting  procedure.  The  coefficients  of  au(R )  were  kept 
constant  at  the  values  derived  from  fitting  the  expression  of 
Eq.  (3)  to  the  scaled  spin-orbit  coupling  constant  (see  Fig.  4). 
We  have  also  performed  fits  with  11  parameters  [8  for  the 
22*  and  2fIM  states,  and  3  for  mu(/?)].  The  optimized  values 
for  au0 ,  au  l,  and  au2  differed  by  less  than  1%  from  the 
scaled  ab  initio  values  and  no  noticeable  improvement  of  the 
quality  of  the  fit  resulted.  The  initial  value  for  the  equilib¬ 
rium  intemuclear  distance  of  the  22*  state  was  set  to  the  ab 
initio  value  of  =3.075  A  derived  from  the  MRCI+0 
(DKH2)  calculation  (see  first  column  of  Table  I).  This  value 
was  supported  by  an  independent  measurement  of  the  rota- 
tionally  resolved  photoelectron  spectrum  of  the  1(1 /2m)  state 
of  Xe£  (Ref.  79)  and  only  changed  slightly  in  the  fit  (see 
Table  VII). 

The  calculated  positions  of  the  vibrational  transitions  of 
the  1(1  /2m),  1(3 /2m),  and  11(1 /2m)  states  using  the  optimal 
parameter  set  summarized  in  Table  VII  are  listed  in  Tables 
VIII,  III,  and  IV  where  they  are  also  compared  with  the 
experimental  values.  The  normalized  root-mean-square  (rms) 
deviation  of  the  best  fit  amounted  to  1.04.  This  number  im¬ 
plies  that  the  quality  of  the  fit  is  excellent.80  The  differences 
between  observed  and  calculated  vibrational  positions  were 
smaller  than  5  cm"1  with  the  exception  of  the  i/+=  1 1  level  of 
the  11(1/ 2m)  state  for  which  the  difference  amounts  to 
5.7  cm"1.  Looking  at  the  behavior  of  the  discrepancies  for 
the  11(1 /2m)  state  (see  Table  IV)  one  recognizes  that  these 
are  much  smaller  for  the  vibrational  transitions  observed  in 
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this  work  ( u+=0-9 )  than  for  those  taken  from  Ref.  32  (u+ 
=  10- 1 1).  A  reanalysis  of  the  data  of  Ref.  32  showed  that  the 
values  for  the  vibrational  levels  i/+=8-ll  of  the  11(1 /2m) 
state  do  not  correspond  to  131Xe132Xe  but  to  another  isoto- 
pomer  of  XeJ,  a  conclusion  which  could  only  be  drawn  with 
the  knowledge  of  the  spectroscopic  information  derived  in 
the  present  study. 

2.  g  states 

Of  the  three  low-lying  gerade  electronic  states  only  the 
1(3/2 g)  state  (48  observed  vibrational  levels)  is  well  charac¬ 
terized  experimentally.  Because  of  the  weakly  bound  nature 
of  the  2X;  state  this  state  is  only  described  by  two  parameters 
(Re  £  and  b ^  ).  The  ab  initio  quantum  chemical  calculations 
suggest  that  the  22*  state  is  repulsive  in  the  intemuclear 
distances  probed  by  this  experiment  (see  Sec.  II).  Therefore, 
and  somewhat  arbitrarily,  we  set  the  value  of  the  equilibrium 
intemuclear  distance  of  this  state  to  the  large  value  of  Re  % 
=  7.5  A  and  kept  it  fixed  during  the  fitting  procedure.  Be¬ 
cause  of  the  lack  of  data  on  the  I(l/2g)  and  11(1  /2g)  states 
the  coefficients  agi(R)  (i=0,l,2)  were  held  constant  at  the 
values  derived  from  fitting  the  expression  of  Eq.  (4)  to  the 
scaled  ab  initio  ag(R)  curve  (see  Fig.  4).  A  total  of  five  pa¬ 
rameters  (1  for  the  2Xg  and  four  for  the  2n^  states)  were 
refined  for  the  g  states. 

The  optimal  parameter  set  for  the  description  of  the  po¬ 
tential  energy  functions  of  the  three  g  states  is  summarized  in 
Table  VII.  The  normalized  rms  deviation  of  the  best  fit 
amounted  to  0.41.  The  low  value  of  this  deviation  must  be 
interpreted  with  some  caution.  Indeed,  the  experimental  data 
set  does  not  contain  enough  information  on  the  1(1/2 g)  and 
11(1/2 g)  states.  The  adjustable  parameter  of  the  2  curve  suf¬ 
fices  to  place  the  position  of  the  only  observed  level  of  the 
11(1/2 g)  state  at  its  observed  position.  The  agreement  be¬ 
tween  the  calculated  and  observed  positions  of  this  level  can 
thus  not  be  considered  as  a  proof  of  the  assignment.  The 
calculated  positions  of  the  vibrational  levels  of  the  1(3/2 g) 
and  11(1/2 g)  states  are  listed  in  Tables  II  and  IV.  A  compari¬ 
son  of  the  calculated  positions  of  the  1(3/2 g)  state  with  those 
from  the  experiment  shows  that  the  maximum  absolute  de¬ 
viation  amounted  to  -1.6  cm-1  for  the  v+=  1 1  level. 

C.  Potential  energy  functions 

The  potential  energy  functions  of  the  three  u  and  the 
three  g  states  of  Xe£  derived  from  the  2X*,  2IIM,  2£*,  and  2FIU 
potential  energy  curves  and  the  /^-dependent  spin-orbit  cou¬ 
pling  matrix  displayed  in  Table  VI  are  presented  in  Figs.  5 
and  6.  The  dissociation  energies  and  equilibrium  intemuclear 
distances  of  all  six  states  are  summarized  and  compared  with 
values  from  other  studies  in  Table  IX. 

1.  Potential  energy  functions  of  the  u  states 

The  potential  energy  functions  of  the  u  states  derived 
from  the  global  fit  are  displayed  as  full  lines  in  Fig.  5.  Com¬ 
pared  to  the  results  of  our  previous  study32  (dashed  lines  in 
Fig.  5),  the  dissociation  energies  of  the  three  u  states  are 
almost  equal  whereas  the  equilibrium  intemuclear  distances 
have  shifted  to  larger  values.  The  increased  bond  lengths  are 


FIG.  5.  Comparison  of  the  potential  energy  functions  of  the  1(1 /2m), 
1(3/2m),  and  11(1 /2m)  states  of  XeJ  derived  in  this  work  (full  line),  in  Ref. 
32  (dashed  line),  and  in  Ref.  31  (circles).  The  dashed  horizontal  lines  indi¬ 
cate  for  each  state  the  positions  of  the  highest  observed  vibrational  level. 

a  result  of  a  new  rotationally  resolved  measurement  of  the 
low-lying  electronic  states  of  XeJ.79  The  potential  energy 
function  of  the  1(1 /2m)  state  is  in  excellent  agreement  with 
that  reported  by  Paidarova  and  Gadea.31  At  intemuclear  dis¬ 
tances  smaller  than  ~4  A  the  potential  energy  curves  of  the 
1(3 /2m)  and  11(1  /2m)  states  derived  here  are  located  higher  in 
energy  than  those  reported  earlier  ’  whereas  the  long-range 
behavior  is  almost  identical.  This  discrepancy  is  likely  to 
originate  from  the  larger  equilibrium  intemuclear  distances 
of  the  22*  and  2IIU  states  and  from  the  /^-dependent  spin- 
orbit  coupling  constant.  The  effect  of  the  /^-dependence  of  a 
is  further  discussed  in  the  Conclusions  section. 

2.  Potential  energy  functions  of  the  g  states 

The  potential  energy  functions  of  the  g  states  derived 
from  the  global  fit  are  displayed  as  full  lines  in  Fig.  6.  The 
most  remarkable  feature  in  this  figure  is  related  to  the  poten¬ 
tial  energy  function  of  the  1(1  !2g)  state  which  is  only  very 
weakly  bound  and  gets  its  character  mainly  from  the  repul¬ 
sive  22*  state.  The  low  dissociation  energy  and  the  large 


FIG.  6.  Comparison  of  the  potential  energy  functions  of  the  1(1  /2g), 
1(3/2 g),  and  11(1  /2g)  states  of  Xe£  derived  in  this  work  (full  line),  in  Ref. 
32  (dashed  line),  and  in  Ref.  31  (circles).  The  dashed  horizontal  line  indi¬ 
cates  for  the  I{3/2g)  state  the  position  of  the  highest  observed  vibrational 
level. 
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TABLE  IX.  Dissociation  energies  De  and  equilibrium  intemuclear  distances  Re  of  the  lowest  electronic  states  of 
Xe£  determined  from  the  potential  energy  curves  displayed  in  Figs.  5  and  6. 


1(1/2  u) 

1(3/2*) 

I(3/2u) 

Kl/2*) 

II(l/2u) 

11(1/2*) 

Reference 

Dt  (cm"1)  7937.4 

7932.7 
7904 
7178 
8691 
8646 
6372 

Re(  A)  3.118 

3.034 
3.114 

3.17 
3.196 

3.18 
3.27 


1829.5 

453.9 

1828.8 

455.7 

1905 

595 

1532 

403 

1163 

527 

1202 

339 

968 

242 

3.683 

4.632 

3.674 

4.390 

3.695 

4.395 

3.19 

4.76 

3.979 

4.382 

3.91 

4.00 

3.97 

4.76 

72.9 

1419.8 

147.8 

1430.2 

162 

1602 

78 

1371 

44 

1690 

Repulsive 

1686 

Repulsive 

968 

6.0% 

4.130 

5.437 

3.985 

5.774 

3.983 

6.35 

4.13 

5.260 

4.059 

Repulsive 

3.97 

Repulsive 

4.23 

436.0 

This  work 

436.1 

32 

368 

31 

242 

42 

307 

13 

Repulsive 

40 

61 

41 

4.199 

This  work 

4.448 

32 

4.773 

31 

5.29 

42 

4.620 

13 

Repulsive 

40 

5.29 

41 

intemuclear  distance  could  explain  why  this  state  has  not 
been  observed  in  spectroscopic  experiments  so  far.  Whereas 
the  dissociation  energy  and  the  equilibrium  intemuclear  dis¬ 
tance  of  the  1(3/2 g)  state  are  comparable  to  those  deter¬ 
mined  in  our  previous  study  the  equilibrium  intemuclear  dis¬ 
tance  of  the  11(1  I2g)  state  has  shifted  to  lower  values. 
Because  of  the  scarce  experimental  data  on  the  g  states,  these 
potential  energy  functions  strongly  rely  on  the  value  of  the 
equilibrium  intemuclear  distance  of  the  2£*  state  obtained  by 
ab  initio  MRCI+0  (DKH2  +  SO)  quantum  chemistry. 

VI.  CONCLUSIONS 

The  potential  energy  functions  of  the  six  low-lying  elec¬ 
tronic  states  of  Xe£  were  derived,  for  the  first  time,  by  ex¬ 
plicitly  considering  the  /^-dependence  of  the  spin-orbit  cou¬ 
pling  constant  as  predicted  by  ab  initio  quantum  chemistry. 
The  effect  of  the  /^-dependence  is  best  manifested  by  com¬ 
paring  the  potential  energy  functions  of  the  1(1  /2m),  1(3/ 2m), 
and  11(1  /2m)  states  calculated  with  an  /^-dependent  spin-orbit 
coupling  constant  au(R )  with  those  calculated  with  an 


/^-independent  spin-orbit  coupling  constant  a=aa tom  using 
the  same  2  and  II  potential  energy  curves  (corresponding  to 
the  parameter  set  for  the  m  states  in  Table  VII).  These  poten¬ 
tial  energy  curves  and  the  corresponding  energy  differences 
are  displayed  in  Fig.  7.  The  influence  of  the  /^-dependence 
gets  significant  at  intemuclear  distances  smaller  than  —5  A, 
where  the  deviation  of  the  value  of  the  spin-orbit-coupling 
constant  from  the  atomic  value  becomes  visible  (see  Fig.  4). 
Compared  to  the  potential  energy  function  of  the  1(3 /2m) 
state  calculated  with  fl=Malom,  that  calculated  with  au(R)  is 
shifted  toward  higher  energies  by  —50  cm"1  in  the  range  of 
intemuclear  distances  probed  by  the  experiment,  whereas 
that  for  the  11(1 /2m)  state  is  shifted  to  lower  energies  by 
60-1 10  cm"1.  The  fact  that  the  /^-dependence  of  a  has  a 
larger  effect  on  the  1(3/2 m)  and  11(1 /2m)  potential  energy 
curves  than  on  the  1(1 /2m)  function  can  be  explained  by  the 
strongly  bound  nature  of  the  curve  and  the  fact  that, 
below  4.5  A,  this  curve  is  separated  from  the  fl  curve  by 
more  than  3000  cm"1. 

The  rms  value  of  1 .04  for  the  best  fit  of  the  u  states  of 
Xe£  suggests  that  the  potential  model  and  the  ^-dependence 


FIG.  7.  Left  panels:  potential  energy  functions  of  the 
(a)  1(1 /2m),  (b)  1(3 /2m),  and  (c)  11(1 /2m)  states  of  XeJ 
calculated  with  /^-dependent  a  JR)  (full  lines)  and 
/^-independent  a=aatom  (dashed  lines)  using  the  same  X 
and  II  potential  energy  curves  (parameter  set  for  the  u 
states  in  Table  VII).  The  dashed  horizontal  lines  indi¬ 
cate  for  each  state  the  positions  of  the  highest  observed 
vibrational  level.  Right  panels:  energy  difference  of  the 
two  curves  displayed  in  the  left  panels  for  the  (d) 
1(1 /2m),  (e)  I(3/2m),  and  (f)  11(1 /2m),  states  of  XeJ. 
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of  the  spin-orbit  coupling  constant  are  adequate  to  reproduce 
the  data  within  the  experimental  uncertainties.  The  experi¬ 
mental  data  on  the  g  states  are  much  less  extensive  and  the 
potential  model  appears  to  contain  too  many  adjustable  pa¬ 
rameters  leading  to  a  rms  value  of  0.41  (see  discussion  in 
Sec.  V).  The  quality  of  the  potential  energy  functions  of  the 
six  low-lying  electronic  states  of  Xe£  could  be  improved  by 
rotationally  resolved  spectroscopic  investigations  of  these 
states  which  would  also  provide  independent  information  on 
the  intemuclear  separations.  In  addition,  measurements  of 
the  two  (1/2  g)  states  would  be  desirable.  At  present,  the 
experimental  data  on  the  g  states  are  too  sparse  for  the  de¬ 
termination  of  reliable  potential  energy  functions. 

The  potential  energy  functions  derived  in  this  work  re¬ 
produce  the  positions  of  all  experimentally  observed  vibra¬ 
tional  levels  within  their  statistical  uncertainties.  This  excel¬ 
lent  agreement,  which  could  not  be  reached  in  earlier 
analyses  of  the  spectra  of  Kr^  (Ref.  30)  and  XeJ  (Ref.  32) 
using  an  /^-independent  spin-orbit  coupling  constant,  is  at¬ 
tributed  to  the  inclusion,  in  the  model,  of  the  /^-dependence 
of  the  spin-orbit  coupling  constant.  The  calculated 
/^-dependence  of  the  g  and  u  functions  can  be  well  described 
by  an  exponentially  decaying  function  [Eq.  (4)]  and  a  Morse- 
type  function  [Eq.  (3)],  respectively.  Preliminary  results  on 
ArJ  and  KrJ  suggest  that  these  functions  are  also  adequate  to 
describe  the  /^-dependence  of  the  spin-orbit  coupling  con¬ 
stant  in  the  other  homonuclear  rare  gas  dimer  ions. 
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